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We study thermal instability and formation of local hot spots in a driven nonlinear NbN su- 
perconducting microwave resonator. White noise injected into the resonator results in transitions 
between the metastable states via a process consisting of two stages. In the first stage, the input 
noise entering the system induces fluctuations in the resonator mode. While, in the second one, 
these mode fluctuations result in phase transitions of the hot spot due to induced temperature fluc- 
tuations. The associated noise-activated escape rate is calculated theoretically, and also measured 
experimentally by means of driving the system into stochastic resonance. A comparison between 
theory and experiment yields a partial agreement. 
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I. INTRODUCTION 

The simple model of noise activated escape of a Brow- 
nian particle over a potential barrier successfully explains 
the basic behavior of a large number of metastable sys- 
tems in nature 1 -. Examples of such systems can be found 
in almost all major fields of science: physics, chemistry, 
biology and even engineering^. For instance, it explains 
biochemical reactions in ac-driven protein 3 , the lifetime 
of zero-voltage state in Josephson junctions 4 *^, the mag- 
netization reversal in nanomagnets£i£, noise-activated 
switching in micro— and nano-^2*ii mechanical oscil- 
lators, and photon-assisted tunneling in semiconductor 
hetrostructures 1 ^. 

A well-known pioneering work on the subject is Kram- 
mer's in 1940. In his seminal papeni^, he derived rela- 
tively simple expressions for the thermally induced es- 
cape rate in a one-dimensional asymmetric double-well 
potential. In general, these escape rate expressions take 
the form of T = T exp(— Ub/ksT), where Ub is the po- 
tential barrier height, ks is Boltzmann's constant, T 
is the temperature (where the limit fc^T <C Ub is as- 
sumed), and Tq is a rate prefactor. Important extensions 
and refinements to this formula aimed either to include 
a wider range of damping regimes 1 ^* 1 ^ 1 ^ or accommo- 
date the solutions to other cases such as nonequilibrium 
systems, have been contributed by many authors over 
the yearsi2^Sii£. Examples of such nonequilibrium sys- 
tems are metastable potentials modulated by determin- 
istic forces^ e. g. the case of stochastic resonanc e) 21 ! 22 , 
or metastable systems subjected to nonwhite nois o 23 ' 24 . 
Moreover, efforts have been invested also in extending 
Krammer's rate theory to describe metastable systems 
in the quantum limiliii 2 ^ 2 ^, where escape is dominated 
by tunneling. 

In the present paper we study the escape rate of 
metastable states of thermally instable superconduct- 
ing stripline resonators both theoretically and experi- 
mentally. In recent studies^Z^ we have experimentally 
demonstrated such instability in NbN superconducting 
resonators. The measured response of the system to a 



monochromatic excitation was successfully accounted for 
by a theoretical model, which attributed the instability to 
a local hot spot in the resonator, switching between the 
superconducting and the normal phases. Nonlinearity, 
according to this model, results due to coupling between 
the equations of motion for both, the mode amplitude 
in the resonator and the temperature of the hot spot. 
The coupling mechanism is based on the dependence of 
both the resonance frequency and the damping rate of 
the resonator on the stripline impedance, which in turn 
depends on the temperature of the hot spot. Moreover, 
we have employed this instability to demonstrate experi- 
mentally intermodulation gain^S, stochastic resonance^, 
self-sustained modulation of a monochromatic driv o 31 i 32 , 
period doubling bifurcation, and noise squeezing 33 . 

In the case of thermally instable superconducting 
stripline resonators, the escape mechanism governing the 
lifetime of the metastable states differs in general from 
many of the examples mentioned above. In this case, the 
input noise induces escape in a two-stage process. The 
direct coupling between the input noise and the driven 
mode leads to fluctuations in the mode amplitude, which 
in turn, induce fluctuations in the heating power applied 
to the hot spot. Consequently, the fluctuating heating 
power, which is characterized by a finite correlation time, 
leads to temperature fluctuations. Escape occurs when 
the temperature approaches the critical value and a phase 
transition takes place in the hot spot. 

The remainder of this paper is organized as follows. 
In Sec. II the steady state solutions of the equation of 
motion for the resonator-mode are derived for the case of 
local heating instability. In Sec. Ill a perturbative ap- 
proach is applied in order to include the effect of thermal 
fluctuations. In Sec. IV an escape rate expression char- 
acterizing the metastable states of the resonator is ob- 
tained. In Sec. V a brief explanation regarding stochas- 
tic resonance measurement is given. While in Sees. V. A 
and V. B stochastic resonance measurement results are 
employed in order to extract some of the transition rate 
parameters characterizing the system. Finally, a brief 
summary concludes this paper in Sec VI. 
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II. STEADY STATE SOLUTIONS 

Consider the case of a superconducting stripline mi- 
crowave resonator weakly coupled to a feedline. Driving 



the resonator by a coherent tone a|" 



injected 



into the feedline, excites a mode in the resonator with an 
amplitude A = Se _JWpt , where uj p is the drive angular 
frequency, b m is a constant complex amplitude propor- 
tional to the drive strength, and B (t) is a complex mode 
amplitude which is assumed to vary slowly on the time 
scale of 1/wp. 



By applying the methods of Gardiner and Collett intro- 
duced in Ref^ one can obtain the following input-output 
relation 



b out = b v. 



(9) 



which relates the output signal a° ut 
off the resonator to the input signal a™ = ^"g-^p* en- 
tering the system. 

Thus, the reflection parameter r in steady state is in 
general given by 



A. Mode Amplitude 

In this approximation, the equation of motion for B 
reads^i 



^ = [i(u; p -oj )- 1 }B- l ^' 1 b m + c m , (1) 

where u Q is the angular resonance frequency, 7 = 71 + 72 , 
where 71 is the coupling factor between the resonator and 
the feedline, and 72 is the damping rate of the mode. The 
term c m represents an input noise with a random phase 



(c m ) - , 
and autocorrelation functions given by 



(c m (t)c ln (t')} = {c m *{t)c m *(t')) = 



(c in (t)c m *{t r )) = Gcj 5{t-t') 



(2) 



(3) 



(4) 



By further assuming a thermal equilibrium condition 
at temperature T e g and a relatively high temperature 
case kBT e s ^> Hujo, one has 



G = 



7 k B T cS 



(5) 



Rewriting Eq. in terms of the dimensionless time 
r = Lu t and using the steady state solution 



i (lu p - w ) - 7 

yields the following compact form 

db , , c m 
— + \b = — , 

CtT LOq 



(G) 



(7) 



where b = B — B^ represents the difference between the 
mode amplitude variable and the steady state solution, 
while A reads 



b out _ 72 - 71 ~ ^ (up ~ ^0) 
b m 72 + 71 - i (cj p - u) ) 



(10) 



which is obtained by substituting B^ of Eq. 10 in the 
input-output relation given by Eq. @ and dividing by 
the input drive amplitude b m . 



by 



B. Heat Balance of Local Heating 



The total power dissipated in the resonator Q t is given 



Qt = ^0272 \B\ 



(11) 



Furthermore, assuming that the resonator nonlinearity is 
dominated by a local hot spot in the stripline resonator, 
and that the hot spot area is sufficiently small in order to 
consider its temperature T to be homogeneous, the heat 
balance equation readsSi 



(12) 



where C is the thermal heat capacity, Q is the power 
heating up the hot spot given by Q = aQt, where 
< a < 1, and W = H (T-T ) is the power of the 
heat transfer to the coolant, which is assumed to be at 
temperature To, where H is the heat transfer coefficient. 

In terms of the dimensionless time r and the dimen- 
sionless temperature given by 



e 



T-Tp 
T c -To 



(13) 



Eq. O reads 



^+5(6-600) =0, 
dr 

where the following quantities have been defined 



(14) 



1 [0J V 



wo) 



(8) 



H 

Cujq 



(15) 
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and 



Mode Fluctuations 



2Ha l2 \B\ 



gC(T c -T ) 

Hence, the steady state solution of Eq. i|14|). reads 



©ooO — 



gC(T c -T ) 



(16) 



(17) 



Moreover, if one further assumes that the fluctuation 
of B around B^, is relatively small, one can rewrite Eq. 
(1141) in the following form 



where 



and / reads 



e-e c 



/ = <?e c 



B n 



(18) 



(19) 



(20) 



In general, when a hot spot is generated or alterna- 
tively diminished in the stripline, it affects the resonator 
parameters ujq, 71, 72, a and may induce as a result 
jumps in the resonance response curve. Moreover, as we 
have already shown in previous publications^*^, most of 
the nonlinear experimental results exhibited by our su- 
perconducting NbN resonators can be modeled to a very 
good extent by assuming a step function dependence of 
the resonator parameters loq, 71, 72, a on the hot spot 
temperature 



u 



72 



U>Q S < 1 

co 0n e > 1 ' 71 



71. © < 1 

7ln © > 1 



72s 
72« 



e < 1 
e > 1 



Q < 1 

e > 1 



(21) 



(22) 



In addition, we have shown that, in general, while dis- 
regarding noise, the coupled equations and l|14|) may 
have up to two different steady state solutions. A su- 
perconducting steady state (S) exists when Oooo < 1, or 
when E < E s , where E s = gC (T c — Tq) /2a s "/2sti- Simi- 
larly, a normal steady state (N) exists when Gooo > 1, or 
when E > E n , where E n = gC (T c — Tq) /2a„72 ra fi. 



III. FLUCTUATIONS 

In this section we assume a nonzero noise term c m (t) 
entering the resonator, thus giving rise to fluctuations 
around the steady state solution. 



In this case the solution of Eq. (0) reads 

b (r) = b (0) e- Xr + — f c m (r') e A ( T '- r )dr' . (23) 
^0 Jo 

For relatively long times jt/loq 3> 1 one gets by using 
Eq. (|2J) a zero mean value of the mode fluctuation b 



(6(r))=0 



(24) 



whereas by using Eqs. J2J, and Q respectively, one ob- 
tains the following autocorrelation functions 



(b( n )b(T 2 )) = (b*(n)b*(T 2 ))=0, (25) 



and 



(b(n)b* (r 2 )) 



- A* |T2 — Ti I 



(26) 



B. Local Heating Fluctuations 

Similarly, the solution of Eq. (|18|1 reads 



where 



(0(r))=0(0)e-° r 
is the mean value of 9 variable and 



Mr) = f T f(r')e^ r '- T )dT' 
Jo 



(27) 



(28) 



(29) 



is the deviation. 

The variance of 8, which is denoted as (A 2 , ( r)), can 
be derived with the use of Eqs. P9"|). J20J, and |26]l. In 

the case of small r, namely the case when gr < 1 and 
|A| t <SC 1, one has to lowest order in r 



<A|(r)> 



3 2 6^ 2Gcj r 2 



\B n 



7 



(30) 



On the other hand, for relatively long times gr 1 one 
finds 



2\ _ G^loo 9"o 



<A^> = 



7 + .9^0 



|-Boo| 2 7 (7 + guof + (u! p - wo) 2 



(31) 
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By taking the square of Eq. (|18(l one obtains 



e+ 9 



d(e 2 

dr 



„2/)2 



where the variable £ is given by 
C(T) = Tr 



(32) 



(33) 



Expressing C,{t) as a sum of a mean value and a devi- 
ation terms in a similar manner to Eq. (|27[) yields 



C(r)=<C(r))+A c (r) . (34) 

To evaluate the variance of C, (r) , which is denoted as 
(^A 2 - (t)Y in the limit of relatively long times we employ 
Eqs. O, lEU, (113, (PI) and get 



Furthermore, by evaluating the integral, substituting in- 
stead of G and 6ooo (given by Eqs. an d fT7|> respec- 
tively), and using the notations 



C 



l (7 + g^o) +(^ P -^o) 2 
2 5^0 (7 + 9^0 ) 



r 



one gets 



2~7r\ 



r = Tq exp 



(uj p -cj ) + 7 (7 + ff^o) 



uj (7 + ff^o) 



(39) 



(40) 



(41) 



where 



2 _ ff 2 G0^ o w o {uj p - ujq) 2 + 7 (7 + 9^o) 



|-Boo| 2 7 (7 + gui f + (uj p - u>of 



(35) 



U 



huso \B a 



(42) 



is the energy stored in the resonator corresponding to the 
steady state amplitude -Boo, and 



IV. ESCAPE RATE 

Escape from S to N states originates from a flux at 
point = 1 (or 9 = 1 — Oqoo) flowing from < 1 to 
> 1, or vise versa for the case of escape from N to S 
states. Thus, the escape rate is given by 



uj I C/(l-0ooo,C)dC 
/o 



(36) 



where / (9, £) is the joint probability distribution func- 
tion of the random variables 9 and £. As was shown 
above, in the limit where jt > 1, the expectation val- 
ues (6*) and (Q vanish. In general, / (9, £) is expected 
to represents a joint normal distribution. Moreover, 6 
and Q become statistically independent as the expecta- 
tion value (9 2 ) becomes time independent. This can be 
readily inferred from the following relation 



(A fl A 



a 



id(e 2 ) 

2 dr 



- (0) (0 



(37) 



Thus, by applying the previous approximations one 
finds 



u>o exp 



r = 



2nJ{A%) A 



(exp - 




U c = hujQ \B C \ 



(43) 



is the mode energy corresponding to the critical ampli- 
tude B c at which ©ooo = 1, namely 



1 — 0ooO — 



2ha l2 \B C \ 2 
gC(T c -T ) 



(44) 



Note that typically in our NbN devices 2 ^ "f/go-io — 10~ 2 . 
Thus, by assuming the limit j/gujQ <C 1, and the reso- 
nance case ojp = ujq, the above expression appearing in 
Eq. (|4*T|l reduces into 



1 = exp 

2tt 



1{Uc-Uqo) 

2 UooksTcff 



(45) 



V. STOCHASTIC RESONANCE 

In order to examine experimentally the escape rate 
expression derived in Eq. Ij41|l . we employed stochas- 
tic resonance technique. Basically, stochastic resonance 
phenomenon demonstrates how a weak periodic signal, 
applied to a nonlinear metastable system, can be ampli- 
fied at the system output with the aid of certain amount 
of zero-mean Gaussian white noise. The amplification of 
the signal occurs when a resonant cooperation is estab- 
lished between the small periodic signal and the white 
noise entering the system. In general, such a coherent 
interaction between the signal and the noise occurs when 
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the angular frequency Q of the signal, which periodically 
modulates the double-well potential of the system, be- 
comes comparable to the escape rate of the metastable 
states in the presence of the white noise. 

The stripline center-layer layout of the NbN resonator 
employed in the measurements is shown at the top-right 
corner of Fig. ^ The resonator was dc-magnetron sput- 
tered on a 34 mmX 30 mmX 1 mm Sapphire substrate 
in an ambient gas mixture of Ar/N2 at room tempera- 
ture. The resonator was patterned using standard optical 
lithography and ion-milling. The resonator thickness was 
set to 2200 A. Additional fabrication process parameters 
are listed in Refill. Whereas modeling and character- 
ization of these nonlinear resonators are elaborated in 
RefA 

The metastable states of the system in our case, are 
manifested by the occurrence of jumps in the resonance 
line shape of the resonator. In Fig. [2] (a) we show 
a reflection parameter measurement of the first reso- 
nance mode of the resonator fa = ujq/2tt ~ 2.575 GHz, 
which exhibits two frequency hysteresis loops forming at 
both sides of the resonant curve as the microwave fre- 
quency is swept in the forward and backward directions. 
In Fig. |21 (b) on the other hand, we show a reflected 
power hysteretic behavior measured at a constant fre- 
quency f p = lo p /2tt = 2.565 GHz which falls within the 
instable region of the fundamental mode, as the input 
power is swept up and down. Thus, in order to drive 
our resonators into metastability, we have applied a co- 
herent microwave signal at frequency f p and input power 
Pq = —21.5 dBm. Moreover, in order to tune the res- 
onator into stochastic resonance condition, we have ap- 
plied a small sinusoidal forcing to the system in the form 
of amplitude modulation, and injected a thermal white 
noise with an adjustable intensity to the resonator port. 

A schematic illustration of the stochastic resonance 
measurement setup used is depicted in Fig. ^ A con- 
tinuous microwave signal at frequency f p is amplitude 
modulated at frequency ffi — H,/2tt = 1 kHz. The modu- 
lated signal which effectively modulates the height of the 
potential barrier is combined with a white noise and fed 
to the superconducting resonator. The reflected signal 
off the resonator on the other hand, is mixed with a lo- 
cal oscillator of frequency f p and measured in the time 
domain using an oscilloscope. Additional information re- 
garding stochastic resonance phenomenon measured in 
these nonlinear superconducting resonators is summa- 
rized in Ref. 30 . 



A. Escape Rate Measurement 

At stochastic resonance condition, the lifetime of the 
metastable states becomes approximately equal to half 
the modulation period. Thus, assuming that the sys- 
tem has two metastable states designated by S u and Sd, 
one obtains at this condition, one metastable state es- 
cape event each half time cycle. This is shown in Fig. 
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FIG. 1: Schematic block diagram of the experimental setup 
used in order to measure stochastic resonance. The microwave 
signal generator and the local oscillator at frequency f p were 
phase locked. The layout of the resonator is shown at the 
top-right corner. 



13 which shows a typical result taken in the time domain 
at stochastic resonance conditions, where the jumps ap- 
pearing in the output signal correspond to alternating 
S u — ► Sd and Sd —> S u transitions. 

The blue dotted line shows the amplitude modulation 
signal, which modulates the escape rates Ti and T2 of 
the transitions Sd — > S u and S u — > Sd respectively. 
Near the minimum (maximum) points of the amplitude 
modulation signal the rate T% ^2) obtains its largest 
value, which is denoted as r m i (r m 2). Let t\ (T2) be the 
difference between the time of the transition Sd — > S u 
(S u —> Sd) and the time at which the corresponding es- 
cape rate obtains its largest value, namely the time at 
which Ti = r m i (r2 = r m 2). The probability density of 
the random variable t\ fa) is denoted by fx fa) (fa fa)). 

An estimate for the escape rates r m i and r m 2 could be 
obtained by measuring the probability densities fa fa) 
and fa fa). As can be seen from Eq. (|A1 2|l in the 
appendix, r m i and r m 2 can be estimated from the ex- 
pectation value and the variance of the corresponding 
random variables T\ and To,. However, a more accurate 
value of the prefactor r m can be obtained by invoking 
Eq. HA3fl an d using the measured probability density 
function / fa) . 

In Fig. 01(a) and (b) we show the measured probability 
densities fi fa) and fa fa) derived from 5000 modula- 
tion cycles sampled in the time domain. The solid line in 
both panels represent a Gaussian function fitted to the 
measured probability density in each case. The transi- 
tion rate Ti ^2) as a function of the random variable t\ 
fa), which is found using Eq. (|A3|) and the Gaussian fit, 
is shown in the inset of Fig. [4] (a) [Fig. f|] (b)]. These 
plots yield also the values T m i ~ 4.6 • 10 5 Hz and F m 2 ~ 
2.7 ■ 10 5 Hz for the transitions Sd — > S u and S u — > Sd 
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FIG. 2: (Color online). (a) Forward and backward fre- 
quency sweeps applied to the first mode of the resonator at 
~ 2.575 GHz . The sweeps exhibit hysteresis loops at both 
sides of the resonance line shape. The plots corresponding 
to different input powers were shifted by a vertical offset for 
clarity, (b) Reflected power hysteresis measured at a con- 
stant angular frequency of lo v — 2ir ■ 2.565 GHz which resides 
within the left-side metastable region of the resonance. For 
both plots the black (dark) line represents a forward sweep 
whereas the green (light) line represents a backward sweep. 



respectively. 



B. Discussion 

In order to obtain theoretical estimates for the escape 
rates r ml and r m 2 corresponding to the Sd — > S u and 
S u — * Sd transitions respectively we rewrite Eq. (|41|l in 
terms of the feedline input power Pi n at the extremum 
of the AM modulation, and the power difference APi n = 
P c — Pi n , where P c is the critical input power, which is 
proportional to U c , 



T m = T exp 



-c 



/AP„ 



V Pin J k B T e f{ 

where, according to Eq. (J), the following holds 



(46) 



Pn 1 - M 



27 



(47) 



Estimates for the model parameters corresponding to 
both Sd — * S u and S u — ► Sd transitions are summa- 
rized in Table 1. Estimates for the coupling parame- 
ter 7 corresponding to the S u and the Sd states have 
been extracted indirectly by fitting Eq. (|10f> to the mea- 
sured r vs. uip curves in the vicinity of the resonance 37 . 



FIG. 3: (Color Online). A typical snapshot of the time do- 
main as the resonator is tuned into stochastic resonance con- 
dition. The solid (green) line represents the reflected mod- 
ulated signal, corresponding to ten modulation cycles out of 
5000 employed in the analysis. The dotted (blue) sinusoidal 
line represents the modulation signal applied to the microwave 
signal generator. 



TABLE I: Calculated and Measured Model Parameters 





Sd > Su 


S u *" Sd 


9 [ID" 3 ] 


1.56 


1.56 


7 [MHz] 


37.6 


18.6 


Pc [dBm] 


-23 


-19.6 


AP m [10~ 6 W] 


0.12 


1.2 


\r\ 


0.55 


0.8 


k B T ctt [fW/Hz] 


1.4 


1.4 


To [Hz] 


8 ■ 10 6 


8.3 • 10 6 


r m [Hz] (calc.) 


7.8 • 10 6 


1.9 ■ 10 6 


r m [Hz] (meas.) 


4.6 • 10 5 


2.7 • 10 s 



Whereas, the cooling parameter g, which is defined in 
Eq. (|15fl . has been estimated using experimentally mea- 
sured material properties of NbN 3 ^* 3 ^^, yielding the 
value g ~ 1.56 • 10 -3 (see Refsi 28 i 32 ). Employing these 
estimates together with the experimental values of P c , 
AP in , r, and substituting in Eq. l|4^|l yield a rough 
estimate for the escape rates r ml ~ 7.8 • 10 6 Hz and 
r m2 ~ 1.9- 10 6 Hz for the S d -> S u and the S u -> S d 
transitions respectively. 

The discrepancy in the values of the escape rates ob- 
tained using the theoretical model as opposed to the ones 
extracted from the experimental data (by about an or- 
der of magnitude) can be attributed most likely to the 
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sition rate is given approximately by 
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FIG. 4: Gaussian probability density functions /i(ti) and 
/2(t"2) fitted to the experimental data which correspond to 
the Sd — > S u transition in panel (a) and to the S u — » Sd 
transition in panel (b). The escape rates Ti and T2 associated 
with both transitions is plotted in the insets of panels (a) and 
(b) respectively as a function of the random time variables ti 
and T2 according to Eq. l|A3jl . 



accumulated errors in the estimated values of the model 
parameters which have been evaluated indirectly. For 
example, the g parameter depends among others on the 
geometry of the hot spot and the thermal properties of 
the deposited NbN film, which are not known precisely. 
Moreover, inaccuracies in the coupling factor 7 may re- 
sult due to some approximations, which were employed 
in the fitting procedure^!. 



T (p) = r m exp 



2 P~Pn 
Pm 



(Al) 



where both r m and k are positive constants. 

The probability distribution function F (r) for a tran- 
sition of the kind Sb — > S a to take place within the time 
interval (— to,r) is given by 



f(t)= r 

J -to 



fit) at 



(A2) 



where / (t) is the corresponding probability density. By 
definition, the following holds 



f(r) 



= r[p(7 



1-F(r) 

The initial condition F (—to) = and Eq. (|A3f> yield 



(A3) 



f(j) = r\p(j)]exp(- J ^ T[p(t)]dt 



(A4) 



Further assume the case where at time t = the func- 
tion p (t) obtains a local minimum p (0) = p m . Near t = 
one has 

p(t) = Pm (l + n 2 t 2 ) +0(t 3 ) . (A5) 
Thus, in the vicinity of t = Eq. IjAip becomes 



r(t) = r m exp(-n 2 n 2 t 2 ) 

and the following holds 



(A6) 



VI. SUMMARY 

In conclusion, a noise-activated escape rate expression 
was derived for the case of a nonlinear superconduct- 
ing resonator having a local-thermal instability. More- 
over, stochastic resonance measurements were exploited 
to experimentally determine the escape rate. A partial 
agreement is found between the theoretical and the ex- 
perimental results. 



/(r) =r m exp {-K 2 n 2 r 2 -V^ 



, 2 2 erf [kTLt) + erf (nflto) 



nil 



(A7) 



Keeping terms up to second order in kQt and assuming 
the case where 



-nflto 



> 1 , 



(A8) 



allow approximating the probability density / (r) by 



APPENDIX A: TRANSITION LIFETIME 

Consider a system which has in general two metastablc 
states designated by S a and Sb and assume that at time 
t = —to the system is in state Sb, where to > 0. The 
transition rate T of the process Sb — > S a depends on an 
externally applied time varying parameter p(t). Further 
assume that for p close to some fixed value p m the tran- 



J{t) = -t= exp 

V7T 



-k 2 ^ 2 I r 



(A9) 



In this approximation the random variable r has a normal 
distribution function with a mean value 



r„ 



2K 2 fl 2 



(A10) 



8 



and a variance 
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